submitted to the Astrophysical Journal 



On Iron Enrichment, Star Formation, and Type la Supernovae in Galaxy 

Clusters 

Michael Loewenstein 1 

Laboratory for High Energy Astrophysics, NASA/GSFC, Code 662, Greenbelt, MB 20771 

loewOmilkyway . gsf c . nasa. gov 

ABSTRACT 

The nature of star formation and Type la supernovae (SNIa) in galaxies in the field 
and in rich galaxy clusters are contrasted by juxtaposing the build-up of heavy metals 
in the universe inferred from observed star formation and supernovae rate histories with 
data on the evolution of Fe abundances in the intracluster medium (ICM). Models for the 
chemical evolution of Fe in these environments are constructed, subject to observational 
constraints, for this purpose. While models with a mean delay for SNIa of 3 Gyr and 
standard initial mass function (IMF) are consistent with observations in the field, cluster 
Fe enrichment immediately tracks a rapid, top-heavy phase of star formation - although 
transport of Fe into the ICM may be more prolonged and star formation likely continues 
to redshifts < 1. The source of this prompt enrichment is Type II supernovae (SNII) 
yielding > 0.1 M Q per explosion (if the SNIa rate normalization is scaled down from its 
value in the field according to the relative number of candidate progenitor stars in the 
3 — 8 Mq range) and /or SNIa explosions with short delay times associated with the rapid 
star formation mode. Star formation is > 3 times more efficient in rich clusters than 
in the field, mitigating the overcooling problem in numerical cluster simulations. Both 
the fraction of baryons cycled through stars, and the fraction of the total present-day 
stellar mass in the form of stellar remnants, are substantially greater in clusters than 
in the field. 

Subject headings: galaxies: clusters: general — galaxies: evolution — stars: formation 
— supernovae: general — X-rays: galaxies: clusters 

1. Introduction 

Rich clusters of galaxies provide a uniquely amenable setting for the study of the complex 
processes and consequences of galaxy formation and evolution. As the largest (> 3 10 14 M ) 
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virialized objects in the universe, with the deepest potential wells, they retain all of their processed 
and unprocessed baryonic matter in the form of stars and (predominantly) ~ 3 — 10 keV plasma 
that are readily studied at infrared/optical/ultraviolet and X-ray wavelengths, respectively. 

Rich clusters are sufficiently vast to be considered representative volumes of the universe and 
are often assumed to be composed of dark matter, stars, and gas in cosmic proportions. However, 
because clusters originate as extreme peaks in the random field of initial cosmological density 
fluctuations, the evolution of their constituent galaxies - and by extension the intergalactic medium 
that they are coupled to - may proceed in a manner that is distinct from their counterparts in the 
field. The resulting biases must be understood, and taken into account, when generalizing from 
cluster-based observables in order to draw conclusions about the universe as a whole. In fact, the 
morphological mix, luminosity function, and star formation history of cluster galaxies all display 
signatures of the effects of the exceptionally dense environments where they form and develop 
(Kuntschner et al. 2002; Croton et al. 2005; Schindler et al. 2005; Romeo, Portinari, & Sommer- 
Larsen 2005). Since galactic outflows associated with the galaxy formation process, in turn, feed 
back on the cluster environment by injecting energy and metals into intergalactic plasma, evidence 
of the cluster /field dichotomy is implanted in the intracluster medium (ICM). 

The distinctions between the evolution of baryons in clusters and in the field may be quantified 
by the juxtaposition of the build-up of heavy metals in the universe inferred from the evolving star 
formation and supernovae rates with the enrichment history of the ICM in rich galaxy clusters - 
the dominant reservoir for metals produced by stars in the ensemble of cluster galaxies. Because, 
under intracluster plasma conditions, abundances of Fe are more easily measured than those of 
other elements they are derived to redshifts > 1 and provide the strongest current constraints on 
the evolution of ICM enrichment. 

Mushotzky & Loewenstein (1997) made an initial investigation along these lines based on their 
analysis of ASCA spectra of clusters out to z = 0.3, and concluded that ICM enrichment was 
inconsistent with concurrent estimates of the star formation rate history derived for field galaxies 
(see, also, Madau, Pozzetti, & Dickinson 1998; Lin k Mohr 2004; Calura & Matteucci 2004). 
Re-examination of this issue is timely given the extension to higher redshift of ICM Fe abundance 
measurements made possible with XMM-Newton (Tozzi et al. 2003; Hashimoto et al. 2004; Maughan 
et al. 2004), and direct estimates of the evolution of supernova rates to comparable redshift in the 
field (Strolger et al. 2004; Dahlen et al. 2004; Barris & Tonry 2006) and in clusters (Gal- Yam, Maoz, 
& Sharon 2002) that constrain the level of supernova metal production per unit star formation. 

Star formation with a standard initial mass function (IMF) produces approximately one core 
collapse supernova (hereafter, SNII) per hundred solar masses of stars formed, while recent empirical 
estimates for the specific rate of Type la supernovae (SNIa) are an order of magnitude lower. For 
typical adopted Fe yields of 0.07 M Q for SNII (slightly higher than the empirical estimate of 
Elmhamdi, Chugai, & Danziger 2003) and 0.7 M for SNIa (Iwamoto et al. 1999), the estimated 
universal fraction of baryons in stars of ~ 0.06 (Fukugita & Peebles 2004) thus implies an Fe 
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enrichment per baryon of ~ 0.065/(1 — R) relative to solar (solar Fe mass fraction = 0.0013) that 
originates in comparable portions from the two classes of supernovae, where R is the stellar mass 
loss fraction integrated over the age of the universe. Since R ~ 40% for a standard IMF (see below), 
one expects a mean universal Fe enrichment per baryon to about one-tenth solar. While this may 
be consistent with the low-redshift (non-cluster) IGM (Prochaska 2004; also, see §4.1 below), it 
falls short by a factor of ~ 4 in the ICM (see, also, Pagel 2002). A simple way to account for this 
shortfall is to invoke a top-heavy IMF (David, Forman, & Jones 1991; Arnaud et al. 1992; Elbaz, 
Arnaud, & Vangioni-Flam 1995; Matteucci & Gibson 1995; Loewenstein &; Mushotzky 1996; Gibson 
& Matteucci 1997) to enhance the formation efficiency of SNII, and perhaps SNIa. Additional 
variations (some linked to the IMF) include the following: (1) a higher average star formation 
rate (more SNII and SNIa), (2) a higher incidence of SNIa per star formed, (3) a higher average 
Fe yield per (Type la and/or Type II) supernova, (4) a significant enrichment by a pregalactic 
stellar population (e.g., Loewenstein 2001). The relative dearth of supernovae in clusters out to 
the highest redshifts for which the rate is presently constrained highlights the cluster abundance 
paradox, and begs the question of whether one or more of the standard assumptions (e.g., an IMF 
that is universal and invariant over time) break down in the universe in general. 

This paper aims to evaluate the plausibility and consequences of possible explanations for ICM 
Fe enrichment, taking into account constraints based on the characteristics of stellar populations in 
cluster galaxies as well as the observed mean evolution of the cluster SNIa rate and Fe abundance; 
and, to examine the resulting implications for our understanding of the star formation history 
(SFH) in the universe and its dependence on environment. 

New data on the histories of star formation, supernovae rates, and elemental abundances mo- 
tivate an approach that goes beyond consideration of the baryon and metal inventories to explicitly 
model their self-consistent evolution. The framework adopted here for constructing such models, 
as presented in §2, is kept simple to restrict the number of possible parameters while allowing 
for consideration of an extensive range in each parameter, and to maintain the transparency of 
the effects and implications of various assumptions and scenarios. Boundary conditions, source 
functions, and parameters are initially designed as appropriate for the universe as a whole; but; 
clear paths for adaptations that might be required for rich galaxy clusters are provided. Models for 
the chemical evolution of Fe in the ICM are presented in §3, and the requirements on the nature 
of star formation and supernovae in cluster galaxies imposed by observations are indicated and 
summarized. §4 includes discussion of the distinction between clusters and the field with respect 
to their intergalactic media and galaxy populations, and more detailed analysis of galactic outflows 
and Fe enrichment in clusters. Implications for the nature of SNIa and its environmental depen- 
dence, comparisons with results of similar investigations, and possible directions for future model 
enhancements are also presented in this section. Conclusions, with an emphasis on the dichotomy 
between clusters and the field, are summarized in §5. I adopt a topologically flat cosmology with 
Hubble constant H Q = 70 km s -1 , and total and baryonic matter densities, relative to critical, 
£l m = 0.3 and Qb = 0.045, respectively (e.g., Fukugita & Peebles 2004). 
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2. Cosmic Chemical Evolution 

2.1. Basic Equations 

I consider the coupled evolution of the globally averaged densities and metal abundances of 
the three following baryonic categories: (1) stars that includes main sequence and evolved stars, 
substellar objects, and stellar remnants; (2) interstellar gas (ISM) defined as the fraction of 
total matter in the gas phase where stars may form at any time in the history of the universe 
(regardless of location); and, (3) intergalactic gas (IGM) defined as the remaining, inert (e.g., 
non-star- forming) fraction of gas-phase matter 1 . Stars and ISM are coupled via star formation and 
mass return, the ISM and IGM via galactic winds (and, potentially, galactic infall - not explicitly 
considered here). Type la and core collapse supernovae directly enrich the ISM. The evolution 
equations (Tinsley 1980; Matteucci & Gibson 1995; Thomas, Greggio, & Bender 1998) in rest- 
frame time for the comoving densities of the three cosmic constituents - p s tars, Pism, Pigm _ and 
for the mass fractions of their "ith" element - /* s tars 2 ) /*ism> /*igm ~ are as follows: 



dpst&rs . 1 1 \ 

— ^ — = PSF - PMR) \\) 

dPlSM . .. , , 

dt = -PSF + PMR - PGW, (2) 

dt 

^.Tstars _ PSF 



dt Pstars 



= PGW, (3) 
(.ASM ~~ .Tstars)' (4) 



df ISM _ PMR / f j f i \ | PsNIa . PSNII /r\ 

~7, — ~ U stars ~ J ISM) + "~ <~ ~ > \°) 

dt piSM PISM PISM 

and 

df \gm _ PGW , t j f i \ / r n 

-J, — - \J ISM ~ / IGM J' l D J 
dt piGM 

where the source terms psp, pmr, and pew are the comoving global mass density rates of star 
formation, mass return, and galactic wind mass transfer from the ISM to the IGM, respectively; 
and, p*gNia an d P*snii are ^ ne density injection rates into the ISM of element "i" from Type la 
and core collapse (mostly Type II) supernovae, respectively. Since the universe is a "closed box," 



Under these definitions, some (small) fraction of the "IGM" may actually consist of hot galactic halo gas, and 
any gas originally in circumgalactic or intergalactic space that subsequently accretes onto galaxies and forms stars is 
considered "ISM." 

2 Note that f x etaiB includes a contribution from metals locked up in remnants that may be significant (Fukugita & 
Peebles 2004). 
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there are five independent variables; /9 s tars + Pism + Pigm = ^bpcrit, where p CI it = 3H^/8ttG is the 
present-day critical density. Also, the total average baryonic abundance is 

,i _ Pstars/* stars + PISm/*ISM + PlGM.TlGM /-n 

/bar" 5— • {') 

"fcPcrit 



2.2. Functions and Parameters 



The goal of this paper is to explore the chemical evolutionary implications of the increasingly 
complete and detailed inventory of matter, and inferred star formation and metal production rates, 
in the universe; and, to evaluate their applicability to the environment of rich galaxy clusters. As 
such, empirical quantities and constraints are directly applied to the fullest possible extent, and 
pegged as points of departure in modeling ICM enrichment. 



2.2.1. The Stellar Initial Mass Function 



The initial mass function of forming stars (IMF) 4>(m) = dN/dm must be specified to determine 
the normalization of the observed star formation rate and total stellar mass, the stellar mass fraction 
recycled to the ISM via mass loss from evolved stars, and the number of supernova explosions per 
unit mass of star formation. Following Kroupa (2001), I adopt a four-part monotonically decreasing 
piecewise power-law IMF, extending from mi = 0.01 M in the substellar regime to m u = 150 M Q 
(Weidner & Kroupa 2004; Figer 2005), with slope {—dlogN / dlogm) ao = 0.3 below the hydrogen- 
burning mass threshold m# = 0.08 M , slope a.\ = 1.3 between mu and ms = 0.5 M , and slopes 
a.2 between mo and mo and 0:3 above mp. 02 = 03 = 2.3 is adopted as a standard; alternatives 
are assumed to have a single-slope above mg {a.2 = 03), or an additional break at mo > 2mo 
where the slope changes from 02 = 2.3 to Q3. 



2.2.2. Mass Return 

The mass density rate of material recycled from stars to the ISM is 

Pmr(*) = / dm4>(m)Am(m)ns F (t - t ms (m)). (8) 

The turn-off mass m to is implicitly given by t ms (m to ) = At where t ms (m) is the main sequence 
lifetime of a star of mass m (Schaller et al. 1992) and At is time since the onset of star formation. 
Am(m) is the mass returned by stars of main sequence mass m, and is assumed independent of 
time. The star formation rate density nsF = Psf/ § m " lu dm$(m)m. 

Am(m) is derived using appropriate remnant masses, m rcm , in the white dwarf, neutron star, 
and black hole regimes (Prantzos, Casse, k, Vangioni-Flam 1993; van den Hoek & Groenewegen 
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1997; Fryer & Kalogera 2001; Woosley, Heger, & Weaver 2002) as follows: 

Am(m) = m — m reia (m); (9) 

m rcm = 0.45 + 0.119m, m < 8; (10) 

m rcm = 1.4, 8 < m < 21; (11) 

m rcm = 1.4 + 0.717(m - 21), 21 < m < 33; (12) 

^rem = 10, 771 > 33, (13) 

where m is in M Q . 

The integrated fraction of mass formed into stars that is returned to the ISM up to the present 
epoch (t now = 13.47 Gyr) is 

R = Jm *°^ m ILJ (14) 

J m; u dm<j)(m)m 

where mto{t n ow) ~ 0.9 M Q is the present epoch main sequence turn-off mass. Under the instanta- 
neous recycling approximation (IRA: t ms (m) — > for all m), pmr(^) = RpSF- The exact treatment 
(equation 8) is adopted for calculating pmr, since star-formation histories that decline on short 
timescales are considered. However, in deriving equations 4 and 5, the rate of metal transfer to the 
ISM via stellar mass loss at time t is estimated as PMR(^)/ i s tars(*) (rather than by evaluating an 
integral similar to that in equation 8). 



2.2.3. Star Formation 

Estimates of /jsf, corrected for dust, are compiled and compared in Hopkins (2004), and 
fitted to analytic functions in Bell (2005) and Strolger et al. (2004). The star formation rate is 
renormalized to assure that 

f n ° W ,, . /,\ Pstars(^now) Pstars (^form) /i r\ 

L, dtKF{t,= — (i^si — • (15> 

where /o s tars(£now) is the cosmic density in stars observed today, and /9 s tars(^form) the density of zero- 
metallicity (Population III) "seed" stars predating the time, tform; when galaxy formation begins 
and the integration of the set of equations (l)-(6) is initiated. 

2.2.4- Core Collapse Supernovae 

The mass density injection rate of element "i" into the ISM by core collapse supernova of 
massive stars ("SNIP) is 

PSNllW = / dm(f>(m)y SN ii(m)nsF{t - t ms (m)), (16) 
J mi 
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where mi < m < ui2 is the range of SNII progenitor masses and ysNii*( m ) is * ne nucleosyn- 
thetic yield of element "i" from a progenitor of mass m. Instantaneous enrichment is an adequate 
approximation for these short-lived stars, reducing equation (16) to 

PSNIl(i) = ^SNII/5SF(t)(y S NII i ), (17) 

where Xsnii = f mi m2 dm(j)(m) / J m mu dm(j)(m)m is the number of SNII per unit mass of star for- 
mation, and (?/snii 1 ) = f mi m2 dm<j){m)ysmi % \m) / fm™" dm(f)(m) the mean SNII yield of the "ith" 
element. 



2.2.5. Type la Supernovae 

In light of the multiplicity of theoretical predictions for the SNIa rate (Barris & Tonry 2006, 
and references therein) I adopt the semi-empirical formalism of Strolger et al. (2004) and Dahlen 
et al. (2004), determining the mass density injection rate of element "i" by SNIa into the ISM from 

PSNIaW = "SNIaMs/SNIa*, (18) 

and ^ 

nskia(i) = vl dt'p SF (t')$(t - 0, (19) 

where v is the number of SNIa progenitor systems per unit mass of star formation, ysNia* is the 
SNIa nucleosynthetic yield of element "i" (assumed to be constant), and ^(td) is the normalized 
delay time distribution function (DTDF) parameterized using the Gaussian distribution found by 
Strolger et al. (2004) and Dahlen et al. (2004) to explain the observed evolution of the SNIa rate, 

$(t d ) = (27Ta 2 )- 1 e- [(td - tc)2/2cT2] . (20) 



2.2.6. Star- Formation- Induced Galactic Wind 

I assume that supernova explosions drive outflow of material from the ISM to the IGM and 
that the mass loss rate per unit volume is proportional to the total supernova rate and the ISM 
density, 

PGw(i) = #Gw(rcSNIa(*) + -F^SNIlPSF W)piSM- (21) 

Metal-rich gas is not preferentially ejected, and SNIa and SNII are assumed to contribute to driv- 
ing outflows with equal efficiency. Given an observationally determined present-day ISM density, 
PiSM^now), setting the galactic wind strength, i^GW is equivalent to assuming a value for the initial 
ISM density, piSM(iform)- I calculate this relationship by integrating equation (2) from t now to tf or m- 
In the absence of infall, pisM^form) must be sufficient to account for all the star formation since 
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ifornn Pstars(^now) — /Ostars(^form) i and is bounded above by the total baryon density. That is, allowed 
values of -Kgw are those that yield 

PlSM,min < PlSM^form) < PlSM.max, (22) 

where pisM,mm = PiSM(t now ) + (p s tars(inow) - Pstars (iform)) corresponds to no wind, and pisM.max = 
^bPcrit — Pstars(iform) to a wind with maximum integrated mass outflow ("maximum wind"). 

Although ram-pressure stripping may be important in determining the abundance gradient, 
it is generally thought to account for a small (though significant) fraction of total ICM metal 
enrichment (e.g., Domainko et al. 2006). However, it may play a role in extending the epoch of 
mass transfer from galaxies to the ICM. I consider an additional galactic outflow term that may 
(at least in part) be associated with this mechanism (see §§3.1.1, 4.2). 

2.2.7. Boundary Conditions 

Equations (1), (2), (4), (5), and (6) are integrated from tf orm to i n ow, with metal-free initial 
conditions: /* stars = /*ism = Pigm = 0. A small seed stellar density is assumed, /O s tars(iform) = 
^iiiPcrit! Pism (iform) is a free parameter within the limits of equation (22). 

2.2.8. Standard and Varying Parameters and Assumptions 

Standard boundary conditions include the formation epoch of the first Population II stars 
(iform = 4.66 10 s yr, corresponding to redshift 10 for the adopted cosmology), the (Population 
III) stellar density fhii/Ocrit = 3 10~ 6 pcrit (Ostriker & Gnedin 1996), and the stellar (p s tars(*now) = 
0.00267p cr it) and ISM (/?isM(inow) = 0.00078p cr i t ) densities at the present epoch (Fukugita & Peebles 
2004). 

The IMF parameters are the slope at high mass, «3, and transition mass, mo, to this slope 
from «2 = 2.3. The choice of IMF, in turn, determines the mass return fraction R (equation 14) 
and renormalization factor for the star formation rate (equation 15), and the number of SNII per 
unit mass of star formation, ifsNii- in the standard model, 03 = 2.3, R = 0.393 (the IRA value is 
0.405), and ifsNll = 0.0104 Mq 1 for m\ = 8 M and ni2 = m u = 150 M Q . The renormalization 
factors for the Bell (2005) and Strolger et al. (2004) star formation rate parameterizations are 0.94 
and 0.62, respectively. I adopt the Bell (2005) SFH parametrization as standard. Since the Strolger 
et al. (2004) rate differs most dramatically in shape at low redshift where a relatively small fraction 
of the integrated star formation occurs, our results are insensitive to this choice (the renormalized 
functions are compared in Figure 2; see below). Adoption of the standard -Ksnii and Psf(^) provides 
consistency with observed SNII rate evolution to redshift 1 (Strolger et al. 2004, Dahlen et al. 2004; 
see, also, Figure 3 below). There is sufficient theoretical uncertainty in the SNII synthesis of Fe 
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that the IMF-averaged SNII Fe yield, (j/snii ), is simply left as a parameter with standard value 
0.07 M . 

Standard SNIa parameters include mean delay time t c = 3 Gyr, dispersion a = 0.2t c , and 
normalization v = 1.5 10~ 3 SNIa M© _1 (see, also, Gal- Yam & Maoz 2004, Greggio 2005) that 
provide the best fit to the observed SNIa rate evolution for the Bell (2005) star formation rate 
parametrization (these slightly differ from their values in Dahlen et al. 2004 due to the different 
adopted star formation rate). I focus on Fe in this paper (thus justifying the neglect of non-explosive 
production in equation 5); a SNIa Fe yield ysNia Fe = 0.7 M is adopted as standard. 

For the above standard parameters and functions, one may vary the galactic wind factor from 
K GW = to K GW max , where $7 6 p crit i^GW max = 3.05 10 3 M Q is such that all of the IGM originates 
in the ISM. 

2.2.9. A Note on Varying the IMF 

The standard model IMF (with slope 2.3 above 0.5 M ), if assumed universal, provides mutual 
consistency among observations of the star formation and SNII rates and build-up of stellar mass 
(Bell 2005; Drory et al. 2005; Gwyn & Hartwick 2005), and with the luminosity density of the 
universe (Baldry & Glazebrook 2003). If the slope is significantly flatter, (1) the higher stellar 
mass loss return fraction (Figure 1) implies that a higher star formation rate than is observed 
would be required to produce the observed amount of stars 3 (Figure 2); and, (2) the implied SNII 
rate would be greater than observed (Figure 3a). Conversely, a steeper high mass IMF slope (i.e., 
ct2 = 2.8; see Kroupa & Weidner 2003) significantly underpredicts the observed SNII rate (Weidner 
& Kroupa 2005) unless the transition at ma is pushed to high mass (Figure 3b). 

2.3. Enrichment in the Standard Model 

The chemical evolution of Fe in the standard model (see entries "1N2.3" in Tables 1 and 2) 
provides predictions to be compared with abundance measurements in stars and interstellar gas 
in field galaxies, and in the intergalactic medium. This serves to indicate possible shortcomings 
in the model and provide a baseline for evaluating what manner of extension or variation might 
be required to explain Fe abundances in rich galaxy clusters. Integrated over a Hubble time, the 
total stellar Fe yield 4 (relative to solar) in the standard model is 1.3 (58/42% from SNIa/SNII) 
corresponding to an enrichment of 0.13 averaged over all baryons (Table 2). The model-predicted 
z = SNII and SNIa rates (Table 2) are consistent with observations (Cappellaro, Evans, & Turatto 



3 Neglecting any effects of a different IMF on the observational estimates of the stellar mass. 
4 Note that this is the yield per star formed and is not renormalized to the present stellar mass. 
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1999). 

The distribution of Fe among stellar, ISM, and IGM components as a function of time depends 
on the strength and time-dependence of galactic outflow. The more numerous SNII dominate the 
galactic wind term (equation 21), but SNIa contribute the majority of Fe enrichment - and do so 
with a significant lag with respect to the accumulation of stellar mass and the ejection of SNII- 
enriched material into the IGM. This results in a relatively recent (z < 1) build-up of Fe in the stars 
(Figure 4a) and ISM (Figure 4b) and a substantial Fe mass fraction contained in these components 
at the present epoch. As a result, little more than half of the total Fe production is injected into 
the IGM for a "maximum" wind (Figure 5). A reapportionment that boosts the IGM Fe fraction 
requires extending the outflow duration, a possible mechanism considered in detail below. 

A comparison, focusing on the z > 2 universe, of the enrichment buildup in these models with 
the inventory of metals derived from observations of damped Lya systems, the Lya forest, and 
Lyman break galaxies is made in §4.1 below. 

3. Modeling Fe Enrichment of the Intracluster Medium 

The formalism of the previous section can be directly adapted to investigate the chemical 
evolution of stars, ISM, and IGM (ICM) in clusters - assuming that they are closed boxes - by 
simply rescaling the physical densities by the mean cluster baryon overdensity factor (i.e., the total 
overdensity multiplied by a cluster baryon "bias" factor - see, e.g., Allen et al. 2004), 5 ~ 10 2 . 
Clearly, the standard model must be adjusted in application to clusters as it predicts an IGM Fe 
abundance > 5 times smaller than measured in the ICM. 

There is strong evidence that star formation ensued more quickly in clusters with respect 
to the field, as one might expect in regions of highest initial overdensity (Schindler et al. 2005; 
Romeo, Portinari, & Sommer-Larsen 2005). However, this in itself does not result in enhanced 
enrichment if the integrated star formation and IMF are unaltered. Estimates for the rich cluster 
mass fraction in stars vary, but are typically twice the estimated universal value of ~ 0.06 (see 
§4.4). Nevertheless, the Fe enrichment shortfall persists if the yield per star formed is that of the 
standard model (Pagel 2002; Nagashima et al. 2005a). Moreover, there are claims of still higher 
stellar mass fractions in galaxy groups (Parker et al. 2006; Eke et al. 2005); yet, group IGM Fe 
abundances are not generally higher than in the ICM (though groups may not be closed boxes). 
Therefore, as elaborated on below, a top-heavy IMF may be required to explain the observed 0.4 
solar ICM Fe abundance. 

While much of the star formation in intermediate and low mass galaxies occurs after the peak 
in the universal SFH at z ~ 1, the most massive (elliptical) galaxies were assembled and their star 
formation completed at earlier epochs (Saracco et al. 2004; Bell 2005; Fontana et al. 2004; McCarthy 
et al. 2004; Hammer et al. 2005; Thomas et al. 2005; Caputi et al. 2005; Jimenez et al. 2005; 
Papovich et al. 2006). Since the fraction of stars formed in elliptical galaxies is higher in clusters 
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and star formation is expected to be accelerated in denser environments (Kodama et al. 2004; 
Thomas et al. 2005), I consider an enhancement in star formation, referred to as the "rapid mode", 
with an exponential time-dependence, pxsF oc e~ t /^ Txtn °' x \ These stars may reside in protogalactic 
fragments, or even in isolation in intracluster space (Lin &: Mohr 2004; Zaritsky, Gonzalez, & 
Zabludoff 2004). Star formation in these extended models is characterized by the dimensionless 
parameter tx, the normalizations of the "normal" (i.e., with the same time-dependence as the 
average field star formation rate; see Figure 2) and rapid star formation modes (expressed, in Table 
1, as the present-epoch stellar baryon fractions: /stars and /xstars, respectively), and the IMF of 
each mode (i.e., mp, Q3, mpx, 03X1 m practice, when considering a nonzero exponential star 
formation contribution, 03 = 2.3 is often adopted). For the models discussed in the next section 
tx = 0.1 is chosen. 



3.1. Accounting for Fe, and its Evolution, in Clusters 

3.1.1. General Considerations 

If the ratio of intergalactic gas to stars in clusters is the same as the universal average, ~ 16, 
severe difficulties, as described in §2, emerge in constructing models that reproduce observed cluster 
Fe abundances. This is illustrated in Figure 6, where the solid curve plots the relationship between 
the ratio of the Fe abundance in galaxies (i.e., a mass-weighted average of stars and ISM) to that in 
the ICM versus the mass fraction of Fe in the ICM, assuming the universal IGM fraction of 0.924. 
A conservative upper limit to the ratio of galactic-to-intergalactic Fe abundances of 5 requires that 
> 70% of the Fe reside in the ICM. However, as indicated above for the standard model, if galactic 
winds are driven by star formation while enrichment is more prolonged due to the delay in SNIa 
explosions, then more than half of the Fe is locked up in galactic stars. This results in a galaxy/ICM 
Fe abundance ratio > 12 (Figure 6). This problem is exacerbated if (as expected) star formation 
is accelerated in clusters relative to the field, but is mollified if star formation is more efficient. A 
galactic-to-intergalactic Fe abundance ratio of 5 implies an ICM Fe mass fraction ~ 55% (~ 45%) 
for a stellar mass fraction that is twice (three times) the universal average (see dotted and dashed 
curves if Figure 6). 

To transport sufficient Fe from galaxies to the ICM requires an outflow that is extended in 
time if the standard treatment of SNIa rates and yields, as previously defined, is adopted. I thus 
modify equation (21) by adding one of the two following terms: (1) pcw(t) - ► PGw(t) + PlSM Awind 
(additional exponential outflow term), or (2) pgw(^) — ► PGw(i) + PiSM,min/ At w j n d (additional con- 
stant outflow term), where /OisM.min is defined following equation (22) and i W md and At w i n d are 
additional parameters (that modify the relationship between i^cw and pisM^form))- 
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3.1.2. Results of Representative Models 

The equations of §2 are solved for an extensive multi-dimensional grid of SFHs (/stars, m Di «3, 
ixstars, rn DX , a 3X ) and outflow (i wind or At wind , K GW or - equivalently - pisM(iform)) parameters. 
Conservatively, further consideration is restricted to models that produce ICM Fe mass fractions of 
40 — 80%. Integrated supernovae Fe yields per baryon of ~ 0.5 — 1 relative to solar are thus implied, 
given the 0.4 solar ICM Fe abundance. The outflow strength and evolution must be adjusted to 
release sufficient Fe into the ICM from galaxies, and affects how quickly the high ICM Fe abundance 
accumulates. 

Of course, there is a large parameter degeneracy among models that predict 0.4 solar ICM Fe 
abundance at z = 0. These are constructed and evaluated in the context of recent data on cluster 
Fe abundance and SNIa rate evolution, not necessarily to narrowly determine these parameters, but 
to reveal the general distinctions between these models and those constructed to explain standard 
universal chemical evolution. We will see that none of the models are fully satisfactory, thus 
motivating exploration of additional variations. 

Details of an illustrative cross section of models that are examined in more detail are displayed 
in Table 1. These are labeled by (1) the present-day total stellar mass fraction relative to the 
universal value (there are equal contributions from normal and rapid modes in all hybrid SFH 
models), and (2) whether the SFH follows the normal (N) or rapid (R) mode functional form, or a 
hybrid (H) of the two. They are further denoted by (3) the high mass IMF slope (03 or a^x), where 
the rapid mode slope is used for the hybrid SFH models (except for model lH0.8Wc where 03 = 1.8, 
the normal mode slope is 2.3 for all H-models; Table 1). The break mass, mo or mox, is usually 
2 Mq and introduces no additional ambiguity otherwise. Finally (4) , the designation "Wx" ( "Wc" ) 
denotes whether an extended exponential (constant) outflow term is included, with the normalized 
timescale displayed in Table 1. The choice of outflow form and timescale, as well as the parameter 
PiSM(iform) (expressed in Table 1 as the initial baryon mass fraction in the ISM, /iSM^form))) likewise 
introduce no ambiguities - i.e. an optimal value is selected for each SFH/IMF combination. The 
models appearing above the demarcation line in Tables 1 and 2 follow the standard SNIa and SNII 
treatments explained in §2, and have tf or m corresponding to z = 10. This subset demonstrates the 
diversity of predicted model outcomes and level of consistency with cluster data, and how these 
vary with changes in the input parameters. 

The present-day characteristics of these models are displayed in Table 2, and the evolution of 
the ICM Fe abundance back to z = 1.5 for these models is illustrated and compared to observations 
in Figures 7ab. Since infall is not considered, ISM abundances are significantly higher than stellar 
abundances (see Figure 4; also, Figure 8b below; and, Pipino et al. 2005). Since one might 
expect higher metallicity ISM to form stars more readily and equalize these, I define the "galactic" 
abundance, (Fe) [, as the mass- weighted average of the ISM and galaxy Fe abundances (see Table 
2); (Fe) 1 serves as an upper limit to the stellar Fe abundance. 

The following general statements can be made. (1) As discussed above, a model with the form 
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of the universal SFH, the universal fraction of baryons in stars, and standard transformation from 
star formation rate to supernova enrichment, falls many factors short of providing sufficient Fe 
(e.g., Model 1N2.3; solid curve in Figure 7a). (2) As illustrated by Model 2H1.05 (dotted curve 
in Figure 7b), basic models where the outflow strictly traces the supernova rate generally predict 
fractions of Fe locked up in stars (and ISM) that are too high to yield consistency with present-day 
ICM Fe abundances (for a reasonable galactic abundance, see discussion above). (3) Flattening 
the IMF whilst retaining the standard model SFH form (e.g., Model lN1.05Wc; dotted curve in 
Figure 7a) produces the observed z = ICM Fe abundance if an extended outflow is assumed; 
however, the buildup of Fe is too recent to be consistent with X-ray spectroscopy of clusters over 
redshifts 0-1 (Tozzi et al. 2003; Hashimoto et al. 2004; Maughan et al. 2004). This is also the case 
in other models with insufficient massive star formation at high redshift, including Model lH0.8Wc 
where the universal value for the baryon fraction in stars originates equally in normal and rapid 
star formation modes (short-dashed curve in Figure 7a). 

The observed magnitude and mild evolution of the ICM Fe abundance is fairly well reproduced 
in a variety of models with a top-heavy IMF, prompt star formation, and extended outflow. These 
include Models lR1.05Wc and 2R1.55Wx where all star formation is in rapid mode (long-dashed 
and dot-long-dashed curves in Figure 7a, respectively), and in a number of hybrid ("2H") models 
that include both star formation modes. 

The "2H" series of models (see Tables 1 and 2, and Figure 7b) are those characterized by 
SFHs producing twice the universal average of the present-day baryon fraction in stars (see §4.4), 
with half originating in normal star formation mode with a standard IMF and half in rapid mode 
with a steeper high mass IMF, that yield z = Fe abundances ~ 0.4 solar. These provide fairly 
good explanations of the observed ICM Fe abundances over z = — 1, provided that the IMF is 
sufficiently top heavy and the outflow is of sufficient magnitude and duration (Models 2H1.05Wx, 
2H1.3Wx, 2H1.55Wx). Earlier ICM enrichment in models with exponential winds, compared to 
that in models with constant winds, results in a flatter decline in Fe abundance with redshift and a 
better match to the observed ICM Fe evolution. The effect of varying the supernova wind strength 
over the full physically allowed range is illustrated in Figures 8ab; Fe abundances span a full range 
A(Fe) ICM ~ 0.1 at each redshift. 

The SFHs for Models 1N2.3, lN1.05Wc, 2H1.3Wx, and 2R1.55Wx are shown in Figure 9 (the 
latter two are representative of the best fits to the observed Fe abundance evolution from Figures 
7ab, the former two are included for comparative purposes). SFHs in dual-mode (hybrid) models 
follow the (scaled) field SFH out to z = 1 since a standard IMF is assumed for the normal mode. 
The star formation rates in the models with a rapid component, of course, greatly exceed the 
field value at high redshift; their SFHs resemble those inferred in massive galaxies (Heavens et al. 
2004; Yee et al. 2005). Similar enhancements are realized in the Type II and Type la supernova 
rates (Figures lOab), although the SNIa explosion delay leads to a turnover in the rate at high 
redshift. The solid curves (Model 1N2.3) match the rates in the field, while the errorbars are 
cluster rates measured by Gal- Yam, Maoz, &: Sharon (2002). The latter are converted from units 
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of "SNU" (supernovae per 100 years per 10 10 solar £>-band luminosities) assuming that clusters 
have twice the overall stellar mass-to-light ratio and twice the stellar fraction as in the field, and 
show surprising consistency with the field rate evolution. Large intermediate redshift SNIa rates 
in clusters are predicted in these models where the abundance and evolution of Fe in the ICM is 
explained via SNII and SNIa from a top-heavy IMF enhancement in high redshift star formation, 
and where the SNIa delay-time characteristics and relative normalization are the same as in the 
field. These rates are marginally inconsistent with the empirical cluster estimates, a characteristic 
shared by all models of the type discussed up to this point that successfully explain the ICM Fe 
abundance and its evolution. 

Additional models with empirically motivated variations in formation epoch and various super- 
nova input parameters are described below the demarcation line in Tables 1 and 2, and discussed 
in the following three subsections. 

3.1.3. Reducing the Mean Supernovae Delay Time 

In order to construct models that better agree with the cluster SNIa data, I consider a re- 
duction, from 3 to 0.5 Gyr, in the mean SNIa delay interval (Maoz & Gal- Yam 2004) that may 
be applied to all SNIa (Model 2H1.05WxStl in Tables 1 and 2). Mannucci, della Valle, & Pana- 
gia (2006) argue for a bimodal population of SNIa that includes one subclass that is associated 
with active star formation and characterized by delay times that are insignificant relative to the 
duration of the star formation epoch, and another subclass associated with quiescent secular evo- 
lution that is characterized by a broad distribution extending to long delay times. I also, therefore, 
consider models where the reduction in mean delay time is exclusively applied to SNIa originating 
in rapid star formation mode (Model 2R1.55WxSt - with only rapid mode star formation, Model 
2H1.05WxSt2). The results are shown in Table 2 and Figures llab. Agreement with observations 
of the ICM Fe abundance evolution is significantly improved, and the z < 1 SNIa rate histories - 
that now more closely trace the SFHs - in the hybrid SFH models are in better accord with the 
observational constraints. 



3.1.4- Reducing the Formation Redshift 

Since the heretofore considered rapid star formation models appear to underpredict the present- 
day cluster SNIa rate, I consider a sequence of rapid SFH models where the initiation of star 
formation is moved forward from z(tf oim ) = 10 to z(tf oim ) = 3. The outcomes of three such models 
(2R1.55WxSz, 2R1.55WxStpz, 2R1.55WxStz; Table 1), distinguished by mean SNIa delay times 
(tcx) that are, respectively, 3, 1.5, and 0.5 Gyr are summarized in Table 2 and Figures 12ab. 
The SNIa rate evolution in these models, optimized to match the observed cluster Fe abundance 
evolution (Figure 12a,) is compared to that in models with earlier formation epoch in Figure 12b. 
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Delaying the onset of star formation improves the match to cluster SNIa rate observations for 
t c x = 0.5 Gyr, but worsens it for t c x = 3 Gyr. 

3.1.5. Naturally Reducing the SNIa Rate 

I also construct models with "naturally" reduced rapid star formation mode specific SNIa rate 
- i.e. the scaling is calculated assuming that a universal fraction of 3 — 8 M0 stars forms SNIa 
progenitor binaries. I consider models characterized by mox = mp = 2 M Q , 03 = 2.3, t c = 3 Gyr, 
z(iform) = 10 or 3, t c x = 3 or 0.5 Gyr, and (ysNii Fe ) = 0.07 or 0.10 M ; and, focus here on an 
illustrative subset of "2H" and "2R" models (i.e., where 12% of present-day cluster baryons are in 
stars formed either in equal proportions in normal and rapid mode, or exclusively in rapid mode) 
with asx = 1-05 and 1.55, respectively. The corresponding rapid mode SNIa rate reduction factors 
are 0.26 and 0.66. 

In order for these models with reduced early SNIa enrichment to match the z > 0.3 ICM Fe 
abundance evolution as well as, e.g., the models presented in Figures 11a and 12a, an increase in 
mean SNII Fe yield to (at least) (ysNii Fe ) = 0.10 M Q is required (Figure 13a). Because the early 
ICM Fe enrichment is now dominated by SNII, the ICM Fe evolution is less sensitive to the choice 
of t cX (for (y S Nii Fc ) = 0.07 M , the fall-off of ICM Fe with redshift in 2H models is steeper if 
t c x = 3 Gyr than if t c x = 0.5 Gyr). Since the differences in the SNIa rate history emerge beyond 
the redshift where observed constraints are available (Figure 13b), the SNIa delay time in clusters 
is unconstrained by observations if the rapid star formation in cluster galaxies was characterized 
by a flat IMF and high SNII Fe yields. Note that with most of the Fe enrichment due to SNII in 
the rapid star formation mode, an extended outflow phase is no longer required - provided that the 
initial ISM fraction, /iSM(iform), is large (e.g., model 2R1.55Stnz). Inclusion of a prolonged outflow 
period reduces the required value of /iSM(iform)- 

3.1.6. Model Result Overview 

A simple model for the chemical evolution of Fe in the universe may be constructed using 
parameters drawn from the observed star formation, and Type II and Type la supernova, rates in 
the field population of galaxies that are mutually consistent for a standard IMF. For Fe yields of 
0.07 M per SNII and 0.7 M per SNIa, the baryons in the universe are enriched to an average 
Fe abundance of 0.13 solar. The universal ratio of galactic-to-intergalactic baryons constrains the 
fraction of Fe locked up in stars and cold gas such that stellar Fe abundances are within a factor of 
two of solar with the precise value dependent on the magnitude and timing of galactic outflows. The 
prolonged epoch of star formation coupled with the SNIa delay shifts much of the Fe enrichment 
to relatively recent (z < 1) epochs, as indicated by recent observations (Lamareille et al. 2006). 

If one naively assumes that rich galaxy clusters are pure representative samples of the universe, 
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this global chemical evolution model directly scales to clusters and Fe abundances in the ICM are 
badly underestimated. However, in order to better understand the differences in the nature of 
star formation in cluster and field galaxies, it is instructive to adopt the field model as a point 
of departure for investigating the variations that are required to reproduce observed z ~ Fe 
abundances of ~ 0.4 solar and mild [z = — 1) ICM Fe abundance evolution. 

Successful variations must include shifts to both higher redshift in the peak of the star formation 
rate, and to flatter high mass slopes in the IMF for star formation at high redshift - implying that 
conversion of gas into stars is, on average, more efficient in rich clusters than in the field. 

Overall Fe enrichment originates in roughly comparable fractions from SNIa and SNII for 
models where the relationship between SNIa and star formation is the same in clusters and the 
field. Due to the long SNIa delay, outflow of Fe from galaxies extends beyond the epoch of rapid 
star formation in order to be incorporated into the ICM; however, if the outflow is too prolonged, 
ICM enrichment is too recent and in conflict with the mild observed Fe evolution (see discussion 
in §4.2 below). The implied SNIa rates at z ~ 1 are somewhat at odds with observational limits 
for clusters, although the conflict is slightly weaker than inferred by Maoz & Gal- Yam (2004) (see 
discussions in §§4.3 and 4.6 below). Unfortunately, there is as yet no data to provide a more 
definitive test by offering a comparison with the very high z = 1 — 2 cluster SNIa rates predicted 
by these models (see Figure 7b; Maoz & Gal- Yam 2004). 

Additional models with SNIa rate histories that differ from those above, and that generally 
improve the match to the z < 1 evolution in ICM Fe inferred from X-ray data, are constructed 
by (1) shortening the mean delay time in the kernel that transforms the star formation into the 
SNIa rate, (2) decreasing the initial star formation redshift in models where all stars form in rapid 
mode, or (3) reducing the number of SNIa explosions per unit mass of stars formed according to 
the number of intermediate mass (3 — 8 M ) stars whilst increasing the average SNII Fe yield to 
0.1 M . The SNIa evolution in models of these type, with both rapid and normal star formation 
modes, most closely reproduce the observed z < 1 supernova rates. Observations of rapid build-up 
of Fe, and relatively low z < 1 star formation and SNIa rates, in clusters are most readily explained 
in models where heavy element enrichment traces star formation. 



4. Discussion, Implications, Predictions 

4.1. The IGM Versus the ICM 

It is valuable to compare, in more detail, observed constraints on the chemical evolution of the 
universe outside of clusters with the standard model of §2, i.e. a model with star formation history, 
stellar mass buildup, and Type II and Type la supernova rate histories as observed in the field, 
and with standard IMF and SNII and SNIa Fe yields. Such an examination yields insights into the 
chemical evolution of the universe, and reveals distinctions from that peculiar to galaxy clusters. 
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Comparisons for the abundance evolution to z = 6 of the stellar, ISM, and IGM components are 
shown in Figure 14. Allowing for an Fe-to-a-element ratio as low as one-third, the IGM abundance 
evolution in models with wind parameter greater than that corresponding to ^iSM(iform) > 0.0076 
(see middle curves in Figure 4 and 5) is consistent with the z = 5 Lya forest lower limit (Songaila 
2001), and with z = 2.5 — 3.5 Lya forest measurements (Schaye et al. 2003; Simcoe, Sargent, Sz 
Rauch 2004) . This implies that > 9% of the z = IGM originates in galaxies and that the average 
stellar Fe abundance at z = 2.5 is ~ 0.05 — 0.1 relative to solar. 

The model Fe abundance averaged over all baryons is also consistent with that measured 
in damped Lya systems at z = 0.4 — 1.5 (argued by Rao et al. 2005 as dominating the mean 
cosmic metallicity in this redshift interval). The model ISM abundances are consistent with the 
z = 2.5 — 3.5 damped Lya measurement if f^isM^form) < 0.028. The standard model (without 
an altered IMF, a temporally extended wind, or adjustment of relative normalizations, delays, or 
yields for supernovae) is in general agreement with all z > 2 non-cluster metallicity observations 
(see, also, Daigne et al. 2005). 

Note that, at least for Fe, there is no evidence of the missing metal problem summarized in, 
e.g., Pettini (2004): the models that assume the empirical z > 2.5 star formation rate are in accord 
with the observations. There are two factors at play here. The primary one is that, although 
the Fe yield per star averaged over a Hubble time is indeed ~ 1.3 solar as generally assumed, 
the average yield prior to z = 2.5 is only ~ 0.6 (see Figure 4) due to the delay in SNIa element 
production. The second factor is that the total ISM density must exceed the value measured in 
damped Lya systems, J7dla(-z = 2.5) ~ 0.0015 (Prochaska & Herbert-Fort 2004). The amount of 
potentially star-forming gas must be sufficient to account for the stars formed over z = — 2.5: 
Oism(2.5) > 0.0027, taking mass return into account and assuming a standard IMF. If one also 
accounts for the mass loss from the ISM necessary to account for the metals observed in the IGM, 
^ism(2-5) > 0.0052 is inferred - more than three times the damped Lya value. Thus, the deficit in 
the detected amount of metals relative to expectations may be better characterized as a "missing 
gas (ISM)", rather than a "missing metal," problem in the z ~ 2.5 universe (see, also, Hopkins, 
Rao, & Turnshek 2005; Prochaska, Herbert-Fort, and Wolfe 2005). 

Conversely, there would be a missing metals problem if the variations considered above needed 
to explain cluster Fe abundances and its evolution were applied to the field, as these were specifically 
introduced to produce a prodigious rapid enrichment. Models where half the stars in the field were 
formed from a flat-IMF, rapid star formation mode with either reduced SNIa delay, or increased 
SNII yields and suppressed SNIa rate, overpredict z = 2.5 Fe abundances outside of clusters by 
> 10. Star formation is fundamentally different in clusters and the field (see §4.4, below). 

A full consideration of limits on departures of the standard model consistent with observations 
outside of clusters is beyond the scope of this work. Given the above conclusions, it is perhaps 
surprising that field and cluster galaxy populations are not more distinct in terms of apparent ages, 
mass-to- light ratios, and abundances/abundance ratios - especially since spheroids dominate the 
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stellar mass in both populations (Bell et al. 2003). A general prediction implied by the successful 
cluster models presented here is of higher Fe abundances in cluster galaxies compared to the field 
(even though a higher percentage of Fe is locked up in galaxies in the field - ~ 60%, compared to 
~ 40% in clusters). This effect is mitigated if metals are preferentially lost in outflows associated 
with the rapid star formation. 

4.2. The Magnitude and Timing of Galactic Outflows in Clusters 

I find that > 25% of the ICM must originate in galaxies (i.e the "ISM") in order for the ICM 
to be enriched in Fe to 0.4 solar. The combination of (1) observations of the large amount of 
intracluster Fe and its mild evolution since z = 1, (2) significant limits on the cluster SNIa rate 
since z = 1, and (3) the assumption that Fe in the ICM originates in galaxies that are known to 
form most of their stars at z > 1 constrains the amount, epoch, and duration of matter and metal 
transport from cluster galaxies to the ICM. Thus observations of the evolution of ICM enrichment 
provide a unique diagnostic of the nature of galactic winds associated with galaxy formation in a 
dense environment. 

Since SNII are more numerous than SNIa at all times for all models 5 (see Figures lOab, noting 
the respective scales), they are primarily responsible for the ejection of material from galaxies. The 
SNII rate (and, hence, SNII Fe enrichment) is proportional to the star formation rate that steeply 
declines with redshift in cluster galaxies. Because the mass loss rate as implemented in equation 
(21) is proportional both to the SNII rate and to the mass in interstellar gas that decreases with time 
as it is consumed by star formation and also lost in galactic winds, it falls even more quickly than 
the star formation rate - and, more quickly than the interstellar medium is enriched. Therefore a 
prompt phase of only modestly enriched mass injection into the ICM ensues; and, an additional, 
more prolonged, outflow (see §3.1.1) must be invoked to transport the necessary Fe into the ICM 
- even if the SNIa delay time is reduced to 0.5 Gyr. The sole exceptions to this requirement are 
models where SNII fully dominate the ICM Fe enrichment that may be assumed to undergo prompt 
and very strong (implying very high initial ISM baryon fractions) outflows, although scenarios with 
longer duration and milder outflows (lower initial ISM baryon fractions) are also feasible (see §3.1.5). 
However, in models where the galactic outflow era is overly prolonged, the resulting delay in ICM 
enrichment leads to an underprediction of the observed ICM Fe abundances at z > 0.3. 

Such an additional extended phase of mass loss may be identified with ram-pressure stripping 
of galaxies rather than galactic winds, since simulations indicate that this is the dominant ICM 
enrichment mechanism at z < 1 (Schindler et al. 2005; Domainko, et al. 2006). In the models 
presented here with extended outflow, ~ 25 — 40% of the enrichment is due to this phase, compared 
to ~ 10 — 25% in the simulations. Other possible avenues of prolonging the enrichment timescale 



5 Except at very recent epochs in rapid-only star formation models, when both rates are low. 
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include efficient winds from a galaxy subpopulation (perhaps late-type or dwarf galaxies) distinct 
from that responsible for the prompt outflow phase, or the suppression at early epochs in the 
conversion efficiency of supernovae energy to the kinetic energy of outflow (e.g., if it were more 
efficiently radiated as a result of a denser protogalactic environment). Some, but not all, these 
mechanisms are particular to the rich cluster environment. 

4.3. Implications for Supernova Physics and Cosmology 

As shown in, e.g., Figure 10b the specific cluster SNIa rate is similar to that in the field 
out to z = 1 - the maximum redshift for which observational constraints are available in either 
environment. The inference, drawn above, that cluster star formation is significantly more efficient 
then implies - assuming universal SNIa Fe yields - that SNIa must occur less frequently per unit 
mass of star formed, and/or predominantly explode at z > 1. 

Figures 15ab revisit selected curves from Figures 10b, lib, 12b, and 13b, focusing on the 
z < 1.3 region, and with the SNIa rates on a linear, rather than logarithmic, scale. The solid lines 
show models with standard SNIa normalizations and DTDFs. These are in marginal agreement 
with the observed cluster rate, and predict an order of magnitude increase just beyond the highest 
observed redshift bin. 

For the hybrid star formation models, decreasing the mean delay time or the number of SNIa 
per star formed according to the IMF improves the match to the observations and results in less 
extreme behavior just beyond z = 1 (Figure 15a), while simultaneously providing models that more 
precisely fit the ICM Fe data (assuming, in the reduced-normalization models, that the SNII Fe 
yield is > 0.1 M Q ). In contrast, the SNIa rate in models where star formation exclusively occurs 
in rapid star formation mode generally decline too steeply at low redshift. 

There is independent justification for characterizing SNIa explosions associated with the rapid 
star formation in cluster galaxies with relatively short delay times. Physically plausible double 
degenerate binary system progenitor models for SNIa with short delay times (Belczynski, Bulik, & 
Ruiter 2005) apparently find realization in regions of active star formation. Evidence includes the 
measurement of higher SNIa rates in blue, relative to red, galaxies (Mannucci et al. 2005); and, 
of higher rates in early- type galaxies when they are radio loud (Delia Valle et al. 2005). This is 
consistent with the observation of systematic differences in the characteristics of individual SNIa 
in elliptical and spiral galaxies (Delia Valle et al. 2005, Garnavich & Gallagher 2005). Finoguenov 
et al. (2002) also suggested that multiple SNIa types are required based on gradients in abundance 
ratios in the ICM. High and persistent Fe abundances, and Fe-to-a ratios, in the nuclei of galaxies 
hosting high redshift quasars also may indicate the presence of short-delay-time SNIa (Dietrich et 
al. 2003, Maiolino et al. 2003). The successful models presented here with reduced SNIa delay 
times in the rapid star formation mode fit in well with these observations, providing a connection 
between star formation in local galaxies and the starbursts responsible for the elliptical galaxy 
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stellar populations that dominate in rich clusters. I have shown that, if the SNIa per star formed 
in such models is determined by renormalizing the field value according to the relative number of 
3 — 8 Mq stars, then the SNII Fe yield must exceed 0.1 M irregardless of the rapid mode mean 
SNIa delay time. However, one need not assume that progenitors of the two proposed subclasses 
of SNIa explode with the same likelihood (Mannucci et al. 2006). 

Even for a given progenitor model, the distribution of delay times is sensitive to the distribution 
of orbital separations and mass ratios, and other factors, that depend on the initial mass function 
and galaxy age and metallicity in a complex and uncertain manner (Greggio 2005). Moreover, 
the evidence for a 3 Gyr mean delay for SNIa in the field has recently been questioned (Barris & 
Tonry 2006). Given these considerations, in combination with uncertainties in the configuration 
of SNIa progenitors and explosion mechanism, one hesitates to suggest that the possibility that 
rapid star formation produce SNIa with shorter delay times casts doubt on the utility of SNIa as 
standard candles (Yungelson & Livio 2000). In any case, since models where SNII dominate Fe 
enrichment are not ruled out, the evolution of intracluster Fe does not constrain the SNIa delay 
time in a model-independent way. Additional consideration of abundance ratios, and particularly 
their evolution, may be more definitive. 

4.4. The Efficiency and Nature of Star Formation in Cluster Galaxies Compared to 

Field Galaxies 

The baryon mass fraction in stars, /b, stars an d star-to-gas ratio, /b, stars //b, gas > m rich galaxy 
clusters may be estimated from several recent observational studies - although uncertainties persist. 
Estimates of the total B-band mass-to-light ratio cluster around 300 M /Lb© with a spread of 
~ 30% (Marinoni & Hudson 2002; Girardi et al. 2002; van den Bosch, Yang, & Mo 2003; Sanderson 
& Ponman 2003), while the stellar mass-to-light ratio for a composite stellar population dominated 
by old stars to the degree appropriate for clusters is ~ 4.5 Mq/Lb© (Marinoni & Hudson 2002; 
note, however, that this assumes a standard IMF). This yields a mass fraction in stars of ~ 0.015 
and /b,stars//b,gas = 1/8 (assuming a cluster baryon bias factor of 0.9) - roughly twice the universal 
ratios. Similar ratios are derived from estimates of the total (Kochanek et al. 2003) and stellar 
(Drory et al. 2004) mass-to-light ratios where luminosities are measured in the K-band (see, also, 
Lin, Mohr, & Stanford 2003). 

The overall implications of these estimates of baryon mass fractions in stars for galaxy clusters 
that significantly exceed the universal value, and their connection to the problem of elemental 
abundances in the ICM, may not be fully appreciated. The modeling approach presented here 
provides a quantitative perspective on some of the distinctions between star formation in the 
field and in clusters that complements studies of stellar populations. In models that successfully 
reproduce observed z = — 1 ICM Fe abundances, a significant fraction of the current mass in 
cluster galaxy stars forms in a rapid mode of star formation characterized by a top-heavy IMF. 
The properties of the populations of Lyman-break and submillimeter galaxies provide independent 
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evidence for such an IMF in the environments of rapid merger-driven star formation at high redshift 
of the kind that dominates in cluster galaxy progenitors (Baugh et al. 2005). Because of the higher 
return fraction in this component, more than two-thirds of the integrated star formation occurs in 
this mode, in accord with the large fraction of cluster stellar mass residing in early-type systems 
(Goto et al. 2003; Postman et al. 2005), and the high formation redshift (Blakeslee et al. 2004; 
Holden et al. 2004; Cimatti et al. 2004; van Dokkum et al. 2004; Chen & Marzke 2004; Thomas et 
al. 2005; Longhetti et al. 2005) and relatively small amount of continuing star formation (J0rgensen 
et al. 2005; Thomas et al. 2005; Romeo, Portinari, & Sommer-Larsen 2005) in cluster galaxies. The 
inferred mass return fraction is ~ 70%, compared to 40% for a standard IMF, implying that star 
formation would be 2-3 times more efficient than would be the case for a standard IMF even 
assuming the same present-epoch stellar mass fraction. If the baryon fraction in stars at z = is 
indeed twice as high in clusters as in the field, and stars form in field galaxies with a predominantly 
standard IMF, then star formation in clusters is 3 - 5 times more efficient. Likewise, ~ 25 — 35% of 
cluster baryons may be inferred to be cycled through stars (this fraction is ~ 10% for the standard 
model). Moreover, although the mass fraction of stars formed that is currently locked up in stellar 
remnants is ~ 13 — 15%, ~ 40 — 60% of the stellar mass today is in that, non-luminous, form 
(compared to ~ 20% for a standard IMF). 

These inferences have implications for constructing and interpreting population synthesis mod- 
els of cluster galaxies, and for estimating their stellar mass-to-light ratios (e.g., Zepf & Silk 1996), as 
well as for evaluating the extent and nature of "overcooling" and feedback in numerical simulations 
of cluster formation (Balogh et al. 2001; Kravtsov, Nagai, & Vikhlinin 2005). It is crucial not to 
neglect mass return in inventories of the integrated star formation, especially in clusters. 

Most, but not all, cluster star formation is required to occur at high redshift. As discussed 
in the previous section, models that provide the best simultaneous fits to the SNIa rate and ICM 
Fe evolution are those with hybrid star formation histories, since the SNIa rate steeply declines at 
z < 1 in models where star formation exclusively originates in rapid mode. This may be taken as 
an indication that the finite (if small) low-redshift cluster SNIa rates imply a non-negligible degree 
of relatively recent star formation activity in cluster galaxies. While the Butcher-Omeler effect 
indicates that recent epoch star formation is suppressed relative to the field by a factor of > 10 
in rich cluster cores, star formation persists outside of the core where late- type galaxies are more 
common (Bai et al. 2006). In any case, since the estimated low- z cluster SNIa rate is only 3a above 
0, it is likely that one could construct hybrid-SFH models with a lower level of recent (z < 0.5 - 
considering the delay in SNIa enrichment) star formation than in those presented here that would 
be equally successful - particularly if one considers a SNIa delay time distribution function that 
falls off less steeply at long delays than the Gaussian function considered here. In fact, if there 
are two subpopulations of SNIa, with DTDFs centered on short (<< 10 9 yr) and long (> 10 9 yr) 
delays, respectively, an exponential distribution for the long-delay SNIa provides a better fit to the 
observed field SNIa rate history (Mannucci et al. 2006). 
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4.5. Origin and Enrichment History of Intracluster Fe 

In the models that, simultaneously, best match the ICM Fe abundance and SNIa rate evolution, 
ICM enrichment traces star formation via either a combination of SNII and short-time-delay SNIa 
(2H1.05WxStl, 2H1.05WxSt2; Figure 11), or domination by SNII (2H1.05WxSn, 2H1.05WxStn; 
Figure 13). For the former, the fraction of stars in SNIa progenitors is assumed to be the same as 
for a standard IMF, even though the total fraction of stars in the 3 — 8 M0 range is lower. Under 
the IRA, the integrated SNIa Fe enrichment per baryon, relative to solar, is 

< F e> te (SMa) . (™) U . ~- Me -0 • <»> 
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where, = /stars + /xstars is the total present-day stellar baryon fraction. The cluster value 



' stars 

inferred from observations is 



(Fe) bar (OBS) = 0.346 (^) (Mp) + 0.135 (^f ) (Fe) gal , (24) 

where /icm is the cluster baryon fraction contained in the ICM; /icm = 0.865 corresponds to 
/stars = 0-12. For (Fe)g a i = 1 — 3 (the successful models here predict values at the higher end of 
this range), (Fe) bar (OBS) « 0.5 - 0.75. 

(Fe) bar (SNIa) is calculated, without assuming the IRA, and plotted in Figures 16ab as a 
function of rapid star formation IMF slope (a^x), assuming vtldx = tub = 2 M , as = 2.3, 
z(tform) = 10, /gt°ars = 0.12, and ysNia Fe = 0.7 M . Evidently, SNIa synthesize a significant fraction 
of cluster Fe only in models characterized by a flat IMF, and with v maintained at its value deduced 
in the field. Moreover, the SNIa Fe yield must be high. While recent multi-dimensional simulations 
(Travaglio et al. 2004; Ropke et al. 2005) predict lower yields, phenomenological studies continue 
to favor higher values of ysNia Fe (Badenes et al. 2006). 

Similarly, the integrated SNII Fe enrichment per baryon, assuming IRA, is 

/ 0.065 \ / j/ S Nii Fc \ (f$g\ ( ggNii \ 
(Fe) bar (SNII) = ^_ j i^^^- j j ^ _ 5 _ j . (25) 

Figures 16ab also include the computation, without assuming the IRA, of (Fe) bar (SNII) for the 
parameters listed above and ysNii Fe = 0.07 M , as well as the summed SN Fe enrichment per 
baryon. In addition, the local IMF slope, and range of (Fe) bar (OBS) estimated above, are indicated. 
These figures emphasize the robustness of the requirement for a flat IMF (unless ysNii Fe > 0.2 M ), 
and illustrate the challenge in producing sufficient Fe in the hybrid star formation models (that 
best match the observed cluster SNIa evolution) if the SNIa progenitor probability for stars in the 
3 — 8 M range is universal. The Fe enrichment per baryon for such models plateaus at ~ 0.55, 
explaining why 7/sNll Fe = 0.1 M (or greater) is favored if v is scaled in this manner (§3.1.5). 

To summarize, the observed frequency of SNIa in clusters (or in the universe as a whole) is 
insufficient to account for cluster Fe enrichment. Therefore, either a separate class of prompt, high 
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Fe yield SNIa that initiate at high redshift explode with an incidence per star formed that exceeds 
what one expects based on the field value and relative number of 3 — 8 Mq stars, or the shortfall is 
made up by SNII. These SNII must primarily originate from rapid star formation with a flat IMF, 
and yield an average of > 0.1 Mq of Fe. If this is the case, SNII dominate cluster enrichment, 
implying a SNII-like pattern of elemental abundances in the ICM and in the oldest cluster galaxies. 
This pattern is fixed at high redshift, unlike its counterpart in field galaxies where Fe enrichment 
is delayed. Enhanced numbers of SNII, and likely gamma-ray bursts as well, in clusters at z > 1.5 
are then expected. 

4.6. Additional Notes on Similar Studies 

Maoz & Gal- Yam (2004) reached similar conclusions with regard to the inadequacy of SNIa, 
and the likely need for a top-heavy IMF, to account for ICM Fe. Their approach was essentially 
as follows. They calculated the number of SNIa required to make up the Fe deficit left after 
subtracting the SNII contribution - assuming a Salpeter IMF (uq = ai = a 2 = a% = 2.35) - from 
what is observed at z = 0. They then distribute those SNIa in redshift according to a particular 
formation epoch tf or m and SNIa delay time distribution - assuming that all stars formed in a burst 
Qbti ^form for comparison with the cluster SNIa rate evolution derived from their observations. In 
agreement with the present work, they find that the specific incidence of SNIa progenitors required 
exceeds the field value by an order of magnitude. Assumptions in their analysis that differ from 
the present one include (1) z(tf OTm ) = 2 or 3, (2) a pure instantaneous burst, (3) a lower stellar 
mass-to- light ratio (implicitly), (4) neglect of the difference between the mass of stars formed over 
a Hubble time and that measured at z = 0. In addition, Maoz & Gal- Yam (2004) adopt a SNIa 
delay time distribution that is both more sharply peaked and with a more prominent tail to long 
delays. Nevertheless, their results are broadly consistent with those presented here in marginally 
allowing scenarios where z(tfo im ) = 3 and the mean delay time is 3 Gyr (Figure 12). The present 
study goes further in explicitly considering a range of IMFs and possible star formation and outflow 
histories to calculate the evolution of ICM enrichment, and in concluding that models with more 
prompt enrichment more accurately predict the flat ICM Fe evolution that is observed. 

Portinari et al. (2004) analyzed the cluster metal inventory at z = 0. An instantaneous burst 
of star formation was assumed, with self-consistent treatment of the stellar mass-to-light ratio and 
mass return, although the stellar Fe abundance was put in by hand. Since chemical evolution is 
not modeled, the observed ICM Fe evolution is not utilized to evaluate scenarios. Models with 
top-heavy IMFs or non-standard SNIa parameters were not explicitly considered. Nevertheless, 
their rejection of a standard IMF is robust and in agreement with the conclusions presented here. 

De Lucia, Kauffmann, h, White (2004) utilized N-body simulations in a ACDM universe, con- 
joined with a parameterized semi-analytic treatment of astrophysical processes in the baryonic 
component such as cooling, star formation, feedback, and galactic winds, to simulate the chemical 
evolution of the IGM and the ICM. Yields appropriate to a top-heavy IMF are adopted, although 
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mass return and the photometric properties of galaxies were computed for a Salpeter IMF. In agree- 
ment with the results of this work, they found that in models tuned to reproduce the observational 
properties of z = galaxies, ICM enrichment was invariably concentrated to high redshifts - a 
direct result of the z = 5 peak in the star formation rate and the assumption of instantaneous 
recycling. 

The semi-analytic models of Nagashima et al. (2005a) included a more detailed treatment of 
chemical evolution, and adopted some of the enhancements listed in §4.7 below, such as a two-phase 
ISM and separate treatment of galaxy disk and bulge components. They fixed the disk and bulge 
IMFs, with the latter characterized by a monotonic flat slope (thus maximizing the SNII/SNIa 
ratio in this population), and also adopt a single prescription for computing the SNIa rate. Neither 
the star formation, nor the SNIa, rate histories were shown or directly compared to observations; 
however, the latter seems to have a mean delay intermediate between those most often considered 
here (i.e. t c = 0.5, 3 Gyr). Yet again, it was concluded that a top-heavy IMF (corresponding to 
a^x = 1) is required to produce sufficient metals - an inference further supported by comparing 
Fe-to-a-element ratios in individual elliptical galaxies with those in galaxy formation simulations 
(Nagashima et al. 2005b). ICM Fe enrichment of the ICM occurs early in their models, as observed, 
with a-element enrichment ensuing even more promptly. 

Contrary to the results of the present investigation, and those described above, Ettori (2005) 
concluded that ICM Fe can be accounted for by enrichment from star formation with a Salpeter 
IMF and SNIa with delay time and normalization as found in the field by Strolger et al. (2004) and 
Dahlen et al. (2004) - similar to model 1N2.3 in the present work that was dismissed as failing to 
produce adequate Fe. Ettori (2005) acknowledged, however, that not enough stars are observed (via 
Lb) for this scenario to be fully coherent. This internal inconsistency can be clarified as follows. 
The supernova parameters adopted in Ettori (2005) result in an Fe yield of 1.05 10 -3 per solar mass 
of star formation compared to the estimated total cluster Fe mass of 8.6 10~ 3 per solar mass of ICM 
-i.e. a mass ratio of stars- formed to gas of 0.82, or a ratio of stars formed to stars observed today 
of 6.3 (for a present-day cluster star-to-gas ratio of 0.12) is implied. The implied return fraction of 
0.84 is grossly inconsistent with the initial assumption of a Salpeter IMF. 



4.7. Model Shortcomings and Future Refinements 

One might be troubled by the high initial ISM baryon fractions, and the high resulting galactic 
Fe abundances 6 , in models that best explain the ICM Fe and SNIa observations. The need for a 
reservoir to provide for the copious ensuing star formation (no additional sources of star-forming 
gas - via, e.g., infall - are subsequently introduced), and to soak up the newly synthesized metals to 



Again, I note that the models do not distinguish between metals locked up in stellar remnants and those obser- 
vationally accessible in still-living stars; the latter can be significantly less metal-rich than the former. 
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be ejected into the IGM drives the requirement for the former. In the context of the simple models 
constructed here, the latter is an inevitable consequence of the early epoch for most star formation 
and the observationally-driven requirement for early enrichment. A number of refinements that may 
prove fruitful in follow-up studies, some of which are relevant to these issues, may be considered. 

In the models presented here, star formation and IGM-enriching outflows originate in the same 
single-component ISM. One may subdivide the ISM into hot and cold phases, with stars exclusively 
forming from the latter and outflows exclusively originating in the former. The two phases can 
exchange mass and be distinctly enriched by supernovae. One may alternatively, or additionally, 
include "bias" parameters to ensure that outflowing and/or star- forming gas is preferentially metal- 
enriched relative to the average in the ISM. These all may affect the relationship between the 
required duration of galactic outflows and the amount of metals locked up in stars. Likewise, 
requirements on the magnitude and timing of outflows may be altered by explicitly considering 
(re)accretion of IGM onto galaxies. Outflows associated with SNII and SNIa may be "decoupled" 
by associating distinct normalizations with each type. 

The number of components may be further expanded in order to, e.g, allow for separate treat- 
ment of bulges and disks within the stellar component as regards star formation rate, IMF, SNIa 
normalization and delay time distribution, outflows/inflows, etc. - rather than the less physically 
motivated division into "field-like" and "cluster-enhancement" components. This additional level 
of detail clarifies the field/cluster dichotomy, and would enable a more detailed evaluation of the 
star formation histories in the models. 

A more sophisticated treatment of supernovae could include yields that depend on progenitor 
mass, thus introducing an IMF-dependence for average yields, or on the metallicity of the progenitor 
population. These are often implemented for SNII, but could be considered for SNIa as well, as 
could a metal dependence of the number of SNIa per star formed. 

A less empirical treatment of star formation could be implemented in a global manner by 
specifying an explicit, parameterized dependence on total ISM mass density constrained by obser- 
vations of field star formation rates and the distinct stellar inventories in clusters and in the field. 
An explicit metallicity-dependent IMF could be considered. The formation and evolution of binary 
star systems may be distinctly and explicitly traced, and granted a role in determining the SNIa 
rate. 

Population III enrichment may be more carefully considered, especially in relation to metal 
build-up at high-redshift. The overall impact of such a component is generally considered minor in 
the IGM out to relatively high redshift (Yoshida, Bromm & Hernquist 2004; Norman, O'Shea, & 
Paschos 2004; but, see Salvaterra & Ferrara 2003), but may be of enhanced significance in clusters 
(Loewenstein 2001). 

Finally, one should consider the abundance evolution of additional elements, such as O, Mg, 
and Si to further constrain the combination of IMF, relative SNII/SNIa normalization, and outflow 
parameters, as well as the average yields themselves. 
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While the above considerations indicate limitations of the present work, the physically and 
observationally motivated, self-consistent models constructed here reveal a number of implications 
and puzzles. By virtue of the simplicity and limited number of free model parameters in this 
approach, one can draw conclusions and identify persistent paradoxes - summarized in the section 
that follows - that are qualitatively model-independent and easily traced to particular phenomenon 
or assumptions. 

5. Summary and Concluding Remarks 

I constructed models for the coupled mass density and Fe abundance evolution of stars, inter- 
stellar (potentially star-forming) gas, and intergalactic (inert) gas in a closed box, with application 
to the universe as a whole and extension to the special environments of rich clusters of galaxies. I 
focused on Fe because of its power as a diagnostic of SNIa astrophysics and the unmatched con- 
straints provided by recent observations of the evolution of the SNIa rate and ICM Fe abundance. 
I adopted empirically based source terms to mitigate the limitations in our understanding of star 
formation and the physics and astrophysics of supernova explosions. I considered the effects of 
varying the star formation rate history, initial mass function, SNII Fe yields, SNIa rate normal- 
ization and distribution of delay times, and strength and duration of galactic outflow on field and 
cluster observables. Confrontation of models with observations of the three components in both 
environments, and over a range in redshift, constrain the mechanisms, and highlight the distinc- 
tions, of metal enrichment in the field and in clusters that depend on the physical characteristics 
of star formation and SNIa explosions. 

5.1. Chemical Evolution of Fe in the Universe 

The average star formation history in the universe is measured out to redshift z ~ 6 (Giavalisco 
et al. 2004), and a SNIa delay time distribution function derived from a comparison with the SNIa 
rate history to z > 1. The results of direct application of these functions - assuming a standard 
IMF, and standard SNIa and SNII Fe yields - to chemical evolutionary models are consistent with 
metallicities measured in damped Lya systems out to z ~ 3.5 and in the IGM to z ~ 5. A standard 
IMF also correctly predicts measurements of the buildup of stellar mass over cosmic time, and 
estimated SNII rates to z ~ 1. A careful accounting of star formation, stellar mass return, and 
galactic winds implies that the ISM density is > 3 times that in damped Lya systems (DLAs) 
at z = 2.5, supporting suggestions that the "missing metal" problem is actually a "missing gas" 
problem - i.e. that there is more potentially star-forming gas than currently measured in DLAs. 
The delay in SNIa enrichment implies that less than 10% of the Fe in field galaxies was in place 
at z = 2 and that the Fe yield per star formed at that redshift was approximately half the current 
value. One expects a negative correlation of the Fe-to-a elemental abundance ratio with stellar 
population age. Although some field star formation may proceed with a flat IMF (van der Wei 
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et al. 2006), drastic departures from the standard model presented here, unless constructed in a 
highly contrived manner, are ruled out by observations and their mutual consistency within the 
standard framework. 

5.2. Chemical Evolution of Fe in Clusters 

The amount of Fe per baryon measured in clusters exceeds that predicted by the successful and 
internally consistent field model described above by more than a factor of four. And yet, surveys fail 
to detect supernovae - in particular SNIa that most efficiently synthesize Fe - in sufficient numbers 
to account for this hyper-enrichment. In this paper I presented, and evaluated the plausibility 
and implications of, possible resolutions to this now well-established (Portinari et al. 2005, and 
references therein) paradox. 

Clusters are typically inferred to have roughly twice the stellar mass fraction as in the field, 
and cluster galaxies are believed to form stars more rapidly. However, these factors are insufficient 
(although the latter is necessary) for reproducing the observed Fe abundance and its mild evolution 
since z ~ 1. Fe must be synthesized in abundance at an early epoch, and must be efficiently 
transported from galaxies to the ICM. Therefore, a significant fraction of star formation must 
proceed rapidly with a top-heavy IMF (see, also, Finoguenov, Burkert, & Bohringer 2003). Star 
formation is 3 — 5 times more efficient in rich clusters than in the field, mitigating the overcooling 
problem in numerical cluster simulations. Both the fraction of baryons cycled through stars, and 
the fraction of the total present-day stellar mass in the form of stellar remnants, are substantially 
greater in clusters than in the field. Assuming that metals are well-mixed in the ISM before ejection, 
~ 40 — 70% of the ICM must originate in galaxies (see, also, Moretti et al. 2003), with the lower 
values generally associated with galactic outflows assumed to continue for longer durations. 

Confirming and extending the conclusions drawn by previous studies (see §4.6), I demonstrated 
that the observed mild amount of evolution in Fe is most accurately reproduced in models where 
the enrichment tracks star formation, implying that synthesis of cluster Fe was dominated by SNII 
(if the SNIa rate normalization is scaled down from its value in the field according to the relative 
number of 3 — 8 M stars) and /or SNIa with short delay times whose progenitors originated during 
a phase of rapid, top-heavy star formation (see, also, Scannapieco & Bildsten 2005). In the former, 
> 0.1 M per SNII is required; in the latter > 10~ 3 M of Fe from SNIa per star formed is required. 

The best matches to the low observed z < 1 cluster SNIa rates are reproduced in models 
where cluster star formation does not exclusively occur at high redshift, but continues to at least 
z ~ 0.5. This is a more prolonged history of star formation than inferred in the central regions 
of rich clusters, and is contingent on the assumption of a steep decline in the SNIa delay time 
distribution function beyond 3 Gyr. 

Observations of high Fe content in the ICM and its persistence to z = 1, in combination 
with constraints on SNIa enrichment from rates measured to similar redshifts, confirms the rapid 
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and efficient nature of star formation in galactic spheroids in clusters, and provides overwhelming 
evidence that this star formation was characterized by a top-heavy IMF that produced Type la less 
efficiently and/or with shorter delay times (on average) than in field galaxies. Rich galaxy clusters 
are, indeed, special environments and should only be treated as "fair samples of the universe" in a 
limited and well-defined context. 

I am grateful to an anonymous referee for carefully considered and constructive comments on 
the original draft of this manuscript. 
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Table 1. Model Input Parameters 
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a Present-day baryon fraction in stars formed in normal mode. 
b Present-day baryon fraction in stars formed in rapid mode. 
c Initial baryon fraction in star-forming gas (ISM). 
d Mean SNIa delay, in Gyr, for normal star formation mode. 
e Mean SNIa delay, in Gyr, for rapid star formation mode. 

f Formation epoch of the first Population II stars moved forward from z = 10 to z = 3. 
g SNIa normalization reduced according to IMF. 
h SNII yield increased from 0.07 to 0.1 M . 
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Table 2. Model Characteristics and Results at z = 
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Note. — Supernova rates, nskn and n-skia are in units of 10 -4 Mpc~ 3 yr~ 1 5~ 1 , where 
5 is the baryon overdensity; /i a (Fe) is the fraction of Fe originating from SNIa, y Fe star s 
the Fe yield per star relative to the solar mass fraction; (Fe) bar the Fe produced per 
total baryon; (Fe) IGM and (Fe) gal the average Fe abundances in the ICM/IGM and in 
galaxies, respectively. 
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Fig. 1. — Integrated return fraction as a function of IMF slope above mo for mo = 0.5, 1, 2, 4, and 
8 M Q (solid, short-dash, long-dash, dot-short-dash, and dot-long-dash line- type, respectively). 
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Fig. 2. — Bell (2005) (solid curves) and Strolger et al. (2004) (broken curves) star-formation rate 
(in M Mpc -3 yr" 1 ) history parameterizations for vtid = 1 M and (top to bottom) ct2 = 0.3, 1.3, 
2.3, and 2.8. The rates are renormalized to yield the measured present-day stellar density when 
integrated over time. (Fukugita & Peebles 2004). 
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Fig. 3. — left(a): Type II supernova rate (in Mpc -3 yr _1 ) evolution corresponding to the curves 
in Figure 2. right(b): Type II supernova rate evolution for Bell (2005) star-formation history 
parameterizations for ct2 = 0.3,1.3,2.3, and 2.8 (solid, dot, short-dash, and long-dash line-type 
respectively). For c*2 = 0.3 and 1.3, the upper (lower) curves denote mp = 1 M Q (mo = 8 M ); 
this is reversed for «2 = 2.8. The histogram shows the observed SNII rates from Dahlen et al. 
(2004) for the redshift intervals z = 0.1 — 0.5 and 0.5 — 0.9 and the field rate from Cappellaro, 
Evans, & Turatto (1999) for z < 0.1 (converted using a mass-to-blue-light ratio of 2.4; Fukugita & 
Peebles 2004); estimated uncertainties are at the ~ 50% level. 
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redshift redshift 



Fig. 4. — left(a): Evolution of stellar Fe abundance for the standard (see text) model. The cases 
with no wind (OiSM^form) = PlSM(*form)/Pcrit = ^stars(inow) + ^iSM^now) = 0.0034), maximum 
wind (^iSM(iform) = — ^ni ~ 0.045) and an intermediate case (r2isM(iform) = 0.0076) are 
denoted by long-dash, dot, and short-dash line-type, respectively. The upper solid curve shows 
the Fe yield per unit star formed (relative to the solar Fe mass fraction), the lower solid curve the 
overall average abundance (or, equivalently, yield per baryon). right(b): Same as (a) for the ISM 
Fe abundance. 
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Fig. 5. — Same as Figure 4 for the IGM Fe abundance (the abundance is for the no- wind case). 
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Fig. 6. — Ratio of galactic Fe abundance to ICM Fe abundance as a function of the present-day 
mass fraction of Fe in the ICM for overall ICM mass fractions 0.924 (solid curve), 0.865 (dotted 
curve), and 0.806 (dashed curve) - corresponding to baryon fractions in stars equal to the universal 
value, twice the universal value, and three times the universal value, respectively (the universal 
value for the ISM baryon fraction is assumed). 
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Fig. 7.— left(a): Evolution of ICM Fe abundance for Models 1N2.3, lN1.05Wc, lH0.8Wc, 
lR1.05Wc, 2H1.3Wx, and 2R1.55Wx (solid, dot, short-dash, long-dash, dot-short-dash, dot-long- 
dash line-type respectively), right(b): Evolution of ICM Fe abundance for Models 2H0.8Wc, 
2H1.05, 2H1.05Wx, 2H1.3Wx, 2H1.3Wc, and 2H1.55Wx (solid, dot, short-dash, long-dash, dot- 
short-dash, dot-long-dash line- type, respectively). Filled circles denote abundances derived from 
ASCA data of kT > 2 keV clusters from Baumgartner et al. (2005) with 90% confidence uncertain- 
ties, open squares those from XMM-Newton data of kT > 5 keV clusters from Tozzi et al. (2003) 
with 2a uncertainties. See S3 and Tables 1 and 2 for model details. 
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Fig. 8. — left(a): Evolution of ICM Fe abundance for Model 2H1.3Wx with variations in i^cw 
ranging from its minimum to maximum allowed value. These correspond to initial ISM baryon 
fractions, /isM(Wn) = 0.30, 0.37, 0.44, 0.51, 0.58, 0.65, 0.72, 0.79, 0.86, 0.93, and 1.0 (lower to 
upper curves). The /iSM(iform) = 0.51 curve is also plotted in Figures 7a and 7b. right(b): Same 
as (a) for galactic abundance (with lower curves corresponding to smaller values of /iSM(iform))- 
The stellar Fe abundance evolution is shown by the broken line. 
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Fig. 9.— Star formation histories for Models 1N2.3, lN1.05Wc, 2H1.3Wx, and 2R1.55Wx (solid, 
dot, short-dash, long-dash line-types, respectively.) Star formation rates are expressed in units 
of <5 _1 Mq Mpc -3 yr _1 where 5 is the baryon overdensity. The Model 1N2.3 curve matches the 
(scaled) observed SFH for field galaxies. 
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Fig. 10. — left(a): Same as Figure 9 for SNII rate histories, with the observed field rate histogram 
reproduced from Figure 3. right (b): Same as (a) for SNIa rate history. Errorbars are observed 
cluster SNIa rates (with la errors) from Table 7 of Gal- Yam, Maoz, & Sharon (2002). Supernova 
rates are in units of <5 _1 Mpc -3 yr -1 . As in Figure 9, solid curves match the (scaled) observed field 
rates. 




Fig. 11. — left(a): Evolution of ICM Fe abundance for models with reduced SNIa delay times 
2H1.05WxStl, 2H1.05WxSt2, 2R1.55WxSt, (solid, dot, short-dash curves respectively). For com 
parison, the evolution for model 2H1.05Wx from Figure 7a is re-plotted (dot-dashed curve) 
right(b): Same as (a) for SNIa rate (in Mpc~ 3 yr^ 1 ) history. 
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Fig. 12. — left(a): Evolution of ICM Fe abundance for models with rapid mode star formation 
initiating at z = 3, rather than z = 10. Shown are models with mean SNIa delay times of 
3 (2R1.55WxSz; solid curve), 1.5 (2R1.55WxStpz; dotted curve), and 0.5 (2R1.55WxStz; short- 
dashed curve) Gyr. right(b): Same as (a) for SNIa rate (in d^ 1 Mpc~ 3 yr _1 ). The z = 10 
counterparts for delay times of 3 (2R1.55Wx; long-dashed curve) and 0.5 (2R1.55WxSt; dot-dash 
curve), Gyr are reproduced. 
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Fig. 13. — left(a): Evolution of ICM Fe abundance with "naturally" reduced SNIa normalization 
and (ysNii Fe ) = 0.1 M Q . In hybrid star formation history models, a%x = 1-05, z(iform) = 10, and 
t c x = 3 (2H1.05WxSn; solid curve) or 0.5 Gyr (2H1.05WxStn; dotted curve); in the rapid star 
formation model, a%x = 1-05, z(t{ OTra ) = 3, and t c x = 0.5 Gyr (2R1.55Stnz; dashed curve). For 
comparison, the evolution for model 2H1.05Wx from Figure 7a is re-plotted (dot-dashed curve). 
right(b): Same as (a) for SNIa rate (in Mpc~ 3 yr" 1 ) history. 
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Fig. 14. — Evolution of the average baryon (solid curve), IGM (dotted curves), stellar (short-dashed 
curves), and ISM (long-dashed curves) standard model Fe abundances for the "maximum wind" 
(IGM: upper curve, stars/ISM: lower curve) and "intermediate" cases shown in Figures 4 and 5. 
The arrow denotes the z = 5 Lya forest lower limit from Songaila (2001); the solid square with 
errorbars the range of z = 2.5 — 3.5 Lya forest measurements from Schaye et al. (2003); Simcoe, 
Sargent, & Rauch (2004), and references therein; the solid circle the z = 2.5 — 3.5 damped Lya 
value with errorbars reflecting the uncertain dust depletion fraction of Fe from Pettini (2004) and 
references therein; and, the open circle with errorbars the z = 0.4 — 1.5 damped Lya systems 
measurement from Rao et al. (2005). The downward broken extensions for the Lya forest limits 
make allowances for the possibility of an a/Fe abundance ratios as high as 3:1. 




Fig. 15.— left (a): SNIa rate (in 1(T 4 5~ l Mpc~ 3 yr^ 1 ) for Models 2H1.05Wx (solid curve), 
2H1.05WxStl (dotted curve), 2H1.05WxSt2 (short-dashed curve), 2H1.05WxSn (long-dashed 
curve), and 2H1.05WxStn (dot-dashed curve), right(b): Same as (a) for Models 2R1.55Wx (solid 
curve), 2R1.55WxSt (dotted curve), 2R1.55WxSz (short-dashed curve), 2R1.55WxStz (long-dashed 
curve), and 2R1.55Stnz (dot-dashed curve). 
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Fig. 16. — left (a): Enrichment per baryon from SNII, SNIa, and their sum as a function of high 
mass IMF slope for hybrid and rapid star formation histories, where the SNIa rate per star formed 
is fixed at its field value. The approximate range of observed values is demarcated by horizontal 
dotted lines. Curves are calculated assuming ysNia Fe = 0.7 M Q and ysNii Fe = 0.07 M Q . right(b): 
Same as (a) with the SNIa rate per star formed scaled according to the IMF. 



